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I. Introduction 

The grid generation problem lends itself to the use of finite elements 
and variational equations. 

(1) Grids are usually generated as smooth solutions to "nice," 
elliptic differential equations--just the equations well suited to 
variational methods. 

(2) The use of smooth finite elements gives the grid a functional 
expression, which can be examined, evaluated, manipulated, and modified 
naturally and cheaply. 

(3) The "grid equations" are chosen for their qualitative character. 
Exactitude of solutions does not matter as long as this is preserved. As 
a result, extremely coarse (cheap) finite elements may generate a grid of 
high quality, if the boundary conditions are well parametrized. 

I succeeded in demonstrating the following; 

(1) Grid-quality solutions of a wide variety of equations--(direct) 
Laplace's, biharmonic, Helmholtz, even nonl inear--can be generated to fit 
reasonable functional boundary conditions in 2D using very coarse rectan- 
gular finite elements, often 6x3 bicubic. I even tried some "wavy" 
operators (with no natural variational expression) to demonstrate the 
method's versatility. I did not try the inverse Laplace equation, but 
I expect no problem but cost. 
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(2) The finite element grids can be refined, locally modified and 
"fine-tuned" using a simple, cheap composition-of- functions approach, 
without having to solve the differential equation repeatedly. 

II. The Finite Elements 

Smooth, rectangular finite elements were used; the ones here are 
bicubic in the interior. For the linear equations, an option of 
cubic boundary conditions with arbitrarily dense nodes was included. In 
the examples shown there are, unless otherwise mentioned, six patches 
"circumferentially" and three "radially," of which only the innermost 
ring of patches "radially" is plotted. 

Where boundary conditions are used, second derivative discontinuity 
in the "circumferential" direction is confined to the ring(s) of elements 
meeting the boundary. 


III. The Variational Expressions 

Both linear and non-linear equations are solved by minimizing a 
variational integral. In the non-linear case there is iteration. 
Equations to fourth order (i.e., expressions squared in the variational 
integral to second order) are treated. In the case of the "wavy" Helm- 
holtz equation H(f) = v^f + k^f = 0, the variation integral for H*H is 
used with Dirichlet ("underdetermined") boundary conditions. 

A variational approach to solving the inverse Laplace's equation is 
known, but was not tried in this research. Other inverse equations (such 
as biharmonic) could also be used. 


Equation 
v^f = 0 
v^f - k^f = 0 
v^f + k^f = 0 
v^f = 0 
v^(f^) = 0 


Variational Integrand 

|vf|2 
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The boundary conditions are always "Dirichlet-like" (Dirichlet for 
second order, Dirichlet or Dirichlet + Neumann for 4th order), which 
allows simple minimization of the integral. 

Since the finite element grids were so coarse, solution of the 
equations was by a one-step symmetric matrix solver. 

IV. Fine-Tuning 

The grids are defined by functional equations x = x(i,j), i and j 
being "counting" variables. But the solutions to the differential 
equations are x = x(s,t). It remains to define s(i,j) and t(i,j). 

The simplest approach is to make them merely linear functions. As 
a matter of fact, I imposed the boundary conditions s(iL,j) = const, 
s(ip,j) = const, t(i,jg) = const, t(i,jj) = const, so that (s,t) lies 
in a box; then I manipulated the interior values through transfinite 
interpolation. If i is "circumferential" the maps are as follows: 


(i) ^ (U) _ 

. v' 

1 = u 

j = c(u,v) 

(unbroken cubic in v) 


TFI 




FE 




Where the F.E. solution has high skewness or is nearly singular, the 
first map allows the requirement on the TFI boundary condition needed to 
mend this to be multiplied by a small constant, avoiding "grid folding." 
A price is paid; the grid comes out looking irregular. It is better to 
avoid the skewness in the FE solution itself, as by using the biharmonic 
equation with perpendicular boundary conditions, 

V. Conclusions 

The finite element approach to grid generation has proved eminently 
successful for linear grid equations. Its coarseness is a very desirable 
characteristic, most encouraging from the point of view of extending it 
to 3D. The wide choice of equations--control led by a small "black box" 
determining the variational integrand--is another asset. So is the 
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method's usability, at Tittle extra cost, with boundary conditions more 
finely specified than the interior finite element grid. 

Iterative solution of nonlinear equations increases the expense by 
more than an order of magnitude, usually requires numerical integration 
of the variational expression (often convenient in the linear cases too), 
and makes it very difficult in general to apply finely-specified boundary 
conditions. I think there are better approaches (see below). 

The "fine-tuning" needs further refinement itself, especially in 
handling variations in the normal velocity of the grid. Choosing the 
grid equation to yield perpendicularity of the finite element grid map 
(possible, for instance, with the biharmonic equation) is a help. Control 
of "circumferential" grid density at chosen locations worked well. 

VI. Future Directions 

Algorithms should be derived to parametrize boundary conditions in 
such a way as to yield good grids using direct (linear) equations, such 
as Laplace's or the biharmonic. The idea is to imitate the Riemann 
mapping, by slowing down where convex (avoiding boundary overlap), and 
speeding up where concave (avoiding "folds" in grid interior). This 
should make use of inverse equations less necessary, and if done right 
should be extendable to 3D. 

The "fine-tuning" algorithms must be refined. They are in principle 
applicable to any grid that can be described as a function. 

For the inverse equations (for instance, Thompson's method), I 
suggest use of linear paneling schemes in (x,y) space. The resulting 
solution can be approximated by (s,t) finite elements simply by solving 
for the (few) internal nodal values of (s,t) and using inverse function 
theorem derivative evaluations. If the solution needs to be iterated, 
linearized variational expressions using this as a starting point should 
be cheap. Such a method might even have uses in linearized flow simula- 
tion, as for cheap streamline tracking. 

The method needs to be extended to 3D, adapted for "block" grids 
(with the equations, if desired, being valid across block boundaries), 
and adapted for vector computers. 
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VII. Figures 

The figures are true representations of the functions they illustrate, 
although some "handwork" was done on some of them to circumvent bugs in 
the evaluation and plotting software. 

All internal finite element grids are six circumferentially by 
three radially, with only the innermost radial layer plotted, unless 
otherwise mentioned. 


Figure 1.- "Heart" - Laplace's equation was approximated, with 
boundary conditions hand-parametrized to give a well- 
conditioned if not unskewed grid. 





Figure 2.- "Joukowski" (Laplace's) - The natural, analytic parametrization 

of a Joukowski airfoil was imitated by a six-node cubic periodic 
spline. Success in avoiding skewness was middling, as the insets 
show. "Radial" velocity at trailing edge was non-zero. 
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Figure 3.- "Joukowski" (Biharmonic) - The biharmonic equation's 


normal derivative condition was used to enforce conformality in 
the limit exactly at boundaries. (This makes the grid Cl near 
the airfoil, and requires the "fine" boundary condition algo- 
rithm.) The insets show its success, and also that "normal" 
velocity at TE is zero. 



Figure 4.- "Helmholtz" - The equation is V^f - k^f = 0, k = .65, with 
0 - s = 36/tt < 6 and 0 - t - 3, radially symmetric boundary condi- 
tions, r(t=0) = 1, r(t=3) = exp(-n). 0 5 t - 1.4 is plotted. 
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Figure 5.- "Wavy Helmholtz" - The equation is V^f + k^f = 0, k = tt/ 3, 
with 0 S s = 36/tt ^ 6 and 0 - t - 6. At t = 0, (x,y) = ( .5+cos9,sin9) , 
and at t = 6, (x,y) = exp(Tr) (cos9,sin0). x includes a wave that 
traverses 1.25 cycle in joining these boundary conditions. The plot 
shows all the grid (except for screen cutoffs), which was solved on 
a 6x6 finite element mesh. 



Figure 6.- Nonlinear - The equatio n is V^( f^) = 0. Thi s 3x3 grid, shown in 
its entirety, approximates x = \/l .75t+l and y = \Js+l , as expected from 
equation and boundary conditions. Slight deviation is visible in x, due 
to the coarseness of the grid and the (2x2) Gaussian integration. 



(A) 



Figure 7.- "Fine tuning," Laplace - (A) = (i,j), 
(B) = (u,v), (C) = (s,t) in the discussion of 
Section IV. The wiggles in (B), due to skew- 
ness, are present only to second order in (C). 



Figure 8.- "Fine tuned" Joukowski (Laplace) - The result of 
composition by the maps of Figure 7 is shown. (The "fine- 
tuning" specifications were deliberately clumsy.) The 
"shock densing" is good and the TE not too bad, but the 
LE shows a nasty glitch where skewness had to be corrected. 
The finite element map is that of Figure 2. 




Figure 9.- "Fine tuning," Biharmonic - Smoother than 
Figure 7, since there is no skewness to be cor- 
rected at the boundary. 



Figure 10.- "Fine tuned" Joukowski (Biharmonic) - Much better 
than Figure 8 at the LE, due to lack of skewness at the 
boundary in the map of Figure 3. The TE is not so good, 
due to adjusting to zero "normal" velocity there, a defect 
more easily correctable than the problem in Figure 8. 




